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BACKGROUND OF THE INVENTION 



1. Field of the Invention 

[0003] The present invention relates to a method for modeling semiconductor 

devices and more particularly to a method for modeling semiconductor devices, such 
as field effect transistors (FET) and high electron mobility transistors (HEMT) for 
relatively accurately determining the physical device characteristics and small-signal 
characteristics to enable the high frequency performance of the device to be forecasted. 

2. Description of the Prior Art 

[0004] HEMT technology provides unparalleled, high-performance 

characteristics at high frequencies (microwave to millimeter wave). As such, HEMTs 
are used in various RF applications. In order to accurately forecast the performance of 
such devices it is necessary accurately model the effect of the components physical 
structure on its high frequency small signal characteristic. Thus, it is necessary to know 
how physical changes to the device will effect device performance in order to determine 
what process changes may be acceptable to improve RF yield product and which may 
be unacceptable which decrease yield. 

[0005] Physical changes in such devices are known to occur as a result of various 

uncontrolled process events, manufacturing equipment changes or intentional process 
enhancement. Currently two methods for modeling the small signal characteristics of 
HEMT devices are known: equivalent circuit modeling; and physical device simulation. 

Equivalent circuit modeling utilizes networks of linear electrical elements to 
model the small signal performance of the device. In the case of HEMT devices, a 
typical equivalent circuit topology is shown in FIG. 1 . This equivalent circuit model 
is known to accurately model measured S-parameters (small signal characteristics) of 
HEMT devices up to 120 GHz. 

[0006] Unfortunately, there is little correlation between the topology of the 

equivalent circuit and the physical structure of the device. The rough correlation of 
each equivalent circuit element to a location and function within a typical HEMT 
structure is shown in FIG. 2. As such, the small signal models are known to perform 
well at modeling measured S-parameters but usually contain model elements that 
drastically diverge from known physical quantities and characteristics. These 



misrepresentations of the physical device are know to manifest as violations of basic 
semiconductor device principles. For example, basic device laws dictate that small 
signal model parameters "scale" in a predictable manner as the periphery is changed. 
However, most models produced by conventional means become less and less accurate 
as scaling is applied. As discussed above, small signal characteristics can also 

be simulated directly from physical device simulators. Such physical device simulators 
utilize comprehensive data about material characteristics and the basic device physics 
to simulate the actual physical location and structure of HEMT devices. Such 
simulators are known to be based upon finite element and Monte Carlo approaches. 
Such analytical tools are adapted to accept input in the form of the device physical 
structure, as generally shown in FIGS. 3, 4 and 5. In particular, these figures show the 
typical cross section and the "epi" stack used for physical simulation of specific device 
structures. In particular, FIG. 3 illustrates a rough scale drawing of a cross section of 
an exemplary HEMT device. FIG. 4 illustrates how the cross section of information 
regarding device structure is input into a known physical device simulator tool, such as 
APDS 1 .0 by Agilent. FIG. 5 illustrates how the epi stack information is input into the 
physical device simulator. FIG. 6 illustrates where the epi stack physically resides 
within the total device structure. 

[0007] Since these tools use a physical structure to simulate performance, the 

correspondence between simulated small signal performance and the device's physical 
characteristics are relatively strong. However the ability of the device simulators to 
accurately model real measured small signal characteristics is relatively inaccurate. An 
example of such results is shown below in Table 1 which provides a comparison of the 
extracted equivalent circuit model elements using equivalent circuit modeling of 
measured S-parameter data and the results as modeled using a physical device 
simulator, for example, using APDS 1.0. 
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Table 1 

Comparison of Modeled Equivalent Circuit Results from Prior Art Methods 



Intrinsic Equivalent 


Equivalent Circuit 


Physical Device 


Circuit Parameter 


Model 


Simulator 


Cgs 


0.227745 pF 


0.1619pF 


•Rgs 


64242 Li 


infinite 12 


Cgd 


0.017019 pF 


0.01019 pF 


Rgd 


133450 Q 


infinite Q 


Cds 


0.047544 pF 


0.044823 pF 


Rds 


160.1791 Q 


736 Q 


Gm 


135.7568 mS 


66 mS 


Ri 


3.034 Q 


2.304 Q 


Tau 


0.443867 pS 


1.433 pS 



[0008] As such, there is a need for a relatively accurate method for relating 

known physical characteristics of a HEMT device to its measured small signal 
characteristics. Specifically accurate methods are needed for producing small signal 
models that are consistent for: measured to model accuracy; physical properties; 
periphery scaling and bias dependence. 

SUMMARY OF THE INVENTION 
[0009] Briefly the present invention relates to a semi-physical device model that 

can represent known physical device characteristics and measured small signal 
characteristics relatively accurately. The semi-physical device model in accordance 
with the present invention uses analytical expressions to model the fundamental electric 
charge and field structure of a HEMT internal structure. These expressions are based 
on the device physics but are in empirical form. In this way, the model is able to 
maintain physical dependency with good fidelity while retaining accurate measured-to- 
modeled small signal characteristics. The model in accordance with the present 
invention provides model elements for a standard small signal equivalent circuit model 
of FET. The model elements are derived from small signal excitation analysis of 
intrinsic charge and electric field as modeled within the device by the semi-physical 
HEMT model. As such, the RF performance can be predicted at arbitrary bias points. 



DESCRIPTION OF THE DRAWINGS 
[0010] These and other advantages of the present invention will be readily 

understood with reference to the following specification and attached drawings 
wherein: 

[0011] FIG. 1 is schematic diagram of an exemplary small signal equivalent 
circuit model for a HEMT device. 

[0012] FIG. 2 is a sectional view of an exemplary HEMT illustrating the rough 
translation of the physical origins for each of the equivalent circuit elements illustrated 
in the small signal circuit model in FIG. 1. 

[0013] FIG. 3 is a cross-sectional view of a HEMT illustrating the regions in the 

HEMT which correspond to the various circuit elements in the small signal equivalent 
circuit model illustrated in FIG. 1 . 

[0014] FIG. 4 is an example of a cross-sectional descriptive input of a physical 

HEMT device structure by a conventional physical device simulation tool. 
[001 5] FIG. 5 is an example of an epi stack descriptive input of a physical HEMT 

device structure for a known physical device simulation tool. 

[0016] FIG. 6 is an example illustrating the location of the epi stack within the 
device structures cross-sectional view. 

[0017] FIG. 7 is an example of a relatively accurate measured-to-model I-V 
characteristics using the semi-physical modeling method in accordance with the present 
invention. 

[0018] FIG. 8 is a elevational view illustrating an epi stack for an exemplary 

HEMT. 

[0019] FIG. 9 is a cross-sectional view of a HEMT for the exemplary epi stack 

illustrated in FIG. 8. 

[0020] FIG. 1 0 is a blown up diagram of the cross-sectional parameters pertaining 

to the T-gate geometry for the exemplary epi stack illustrated in FIG. 8. 
[0021] FIG. 1 1 is a diagram of an electric conductance model used in the semi- 

physical example. 

[0022] FIG. 12 is a Smith chart illustrating the measured vs modeled S- 

parameters SI 1, S 12 and S22 simulated in accordance with the method in accordance 
with the present invention. 



[0023] FIG. 13 is similar to FIG. 12 and illustrates the measured vs. modeled 

values for the S21 parameter. 

[0024J FIG. 14 is similar to FIG. 12 but for the S12 S-parameter. 

[0025] FIG. 15 is a graphical illustration of the semi-physically modeled vs 

measured small signal Gm. 

[0026] FIG. 16 is a graphical illustration of the semi-physically simulated bias 
dependence of the small-signal output conductance Rds. 

[0027] FIG. 17 is a graphical illustration of the semi-physically simulated bias 

dependence of the small signal gate-source and gate-drain capacitance Cgs and Cgd. 
[0028] FIG. 18 is a graphical illustration of the semi-physically simulated bias- 

dependence of the small signal gate source charging resistance Ri. 
[0029] FIG. 1 9 is a graphical illustration of the semi-physically bias dependence 

of the small signal source and drain resistance Rs and Rd. 

[0030] FIG. 20 is a graphical illustration of the measured vs modeled bias 

dependent gain at 23.5 Ghz for a K-band MMIC amplifier. 

[0031] FIG. 21 A and 21B are graphical illustrations of the extracted parameters 

from measured device I-V's for process control monitor testing. 
[0032] FIG. 22 is a graphical illustration of the measured vs semi-physically 

simulated process variation for Gmpk and Vgspk. 

[0033] FIG. 23 is a graphical illustration of the measured vs semi-physically 

simulated process variation for Idpk and Gmpk. 

[0034] FIG. 24 is a graphical illustration of the measured vs semi-physically 

simulated process variation for Imax and Vpo. 

[0035] FIG. 25 is a graphical illustration of the measured/extracted vs semi- 

physically simulated process variation for the small signal equivalent model Rds and 
Gm. 

[0036] FIG. 26 is a graphical illustration of the measured/extracted vs semi- 

physically simulated process variation for the small signal equivalent model Cgs and 
Gm. 

[0037] FIG. 27 is a graphical illustration of the measured vs semi-physically 

simulated physical dependence for Imax as a function of physical gate length. 
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[0038] FIG. 28 is a graphical illustration of the measured/extracted model vs 

semi-physically simulated physical dependence for Rds as a function of physical recess 
undercut width. 

DETAILED DESCRIPTION 
[0039] The present invention relates to a semi-physical device model which can 

be used to simulate RF performance through physically-based device models. The 
semi-physical model is an analytical model based upon empirical expressions that 
model the physics of HEMT operation, hence the terminology "semi-physical". The 
model incorporates real process parameters, such as gate length recess, etch depth, 
recess undercut dimensions, passivation nitrite thickness, and the like. By using 
empirical expressions, the semi-physical model is able to maintain relatively good 
measured to modeled accuracy while accounting for the effects of process variations on 
the device performance. 

[0040] The semi-physical model in accordance with the present invention 

provides model elements for the standard small signal equivalent circuit model or FET 
as illustrated in FIG. 1 . However, unlike conventional methods, the model elements are 
derived from small signal excitation analysis of the intrinsic charge and electric fields 
within the device. As such, the simulated small signal model elements represent a 
relatively accurate physical equivalent circuit description of a physical FET. 
[0041] The general methodology for the semi-physical modeling of intrinsic 

charge, electric conductance and electrical field is as set forth below. First, the 
relationships between the conduction band offsets and electrical permitivities and 
material composition for the various materials in the epi stack are determined. These 
relationships can be performed analytically or by fitting simulated data from physical 
simulators. Subsequently, the basic electron transport characteristics in any of the 
applicable bulk materials in the epi stack are determined. Once the electron transport 
characteristics are determined, the undepleted linear channel mobility is determined 
either through material characterization or physical simulation. Subsequently, the 
Schottky barrier height value or expressions are determined. Once the Schottky barrier 
height value is determined, the semi-physical equations are constructed modeling the 
following characteristics: 
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[0042] Fundamental-charge control physics for sheet charge in the active channel 

as controlled by the gate terminal voltage. 

[0043] Average centroid position of the sheet charge within the active channel 

width. 

[0044] Position of charge partitioning boundaries as a function of gate, drain and 

source terminal voltages. 

[0045] Bias dependence of linear channel mobility and surface depleted regions. 

[0046] Bias dependence of the velocity saturating electric field in the channel. 

[0047] Saturated electron velocity. 

[0048] Electrical conductance within the linear region of the channel, under the 

gate. 

[0049] Electrical conductance within the source and drain access regions. 

[0050] Once the semi-physical equations are determined, the empirical terms of 

the semi-physical modeling equations are adjusted to fit the model I-V (current voltage) 
characteristics against measured values. Subsequently, the empirical terms are 
interactively readjusted to achieve a simultaneous fit of measured C-V (capacitance- 
voltage) and I-V characteristics. Lastly, the empirical modeling terms are fixed for 
future use. 

[0051] By constructing a comprehensive set of semi -physical equations that cover 

all of the physical phenomenon as mentioned above, the physical operating mechanisms 
within a HEMT device can be relatively accurately determined. FIG. 7 illustrates a set 
of relatively accurate measured-to-modeled I-V characteristics for a HEMT using the 
semi-physical modeling discussed herein. In particular, FIG. 7 illustrates the drain-to- 
source current 1^ as a function of the drain-to-source voltage V^ for various gate biases, 
for example, from 0.4V to -1 .0V. As shown in FIG. 7, solid lines are used to represent 
the semi-physical model while the Xs are used to represent measured values. As shown 
in FIG. 7, a close relationship exists between the measured values and the modeled 
parameters. 

[0052] An example of semi-physical modeling for physical device operation in 

accordance with the present invention is provided below. The example utilizes an 
exemplary device as illustrated in FIGS. 8 and 9. Table 2 represents exemplary values 
for the physical cross-section dimension parameters in the model. FIG. 10 relates to a 



blown up T-gate characteristic which is correlated to the parameters identified in Table 
2. 

Table 2 

Values for the Physical Dimension Parameters Input into Device Cross Section 



Layout Parameter 


Units 


Value 


f^atp 1 f»nnth 


La 


fMmi 


0.150 


wing Length 


Lgw 


FMml 


0.520 


Gate Mushroom Crown Length 


Lgmci 


|Mmj 


li Oftft 

u.zuu 


Total Gate Height 


Hg 


[Mm] 


U.ODU 


Gate Stem Height 


Hgstem 


[Mm] 


U.oUU 


Gate Sag Height 


Hgsag 


[Mm] 


0.100 


Gate Cross-Sectional Area 


Gate Area 


[Mm 2 ] 


0.187 


Max Cross-Sectional Area 


MaxArea 


[Mm 2 ] 


0.364 


Total Gate Periph 


Wg 


[Mm] 


200.000 


# Fingers 


N 


[] 


4.000 


trrip-rirain Snapinn 


Dsd 


[Mm] 


1.800 


Gate-Source Spacing 


Dsg 


fMml 


0.700 


Gate-Drain Spacing 


Dgd 


riiml 


i inn 


Gate-Source Recess 


RECsg 


[urn] 


0.160 


Gate-Drain Recess 


RECgd 


[Mm] 


0.240 


Recess Etch Depth 


Hrec 


[A] 


780.000 


SiN Thickness 


Hsin 


[A] 


750.000 


Gatefeed-Mesa Spacing 


Dgfm 


[Mm] 


2.000 


Gateend-Mesa Overlap 


Dgem 


[Mm] 


2.000 


Finger-Finger Spacing Thru Drain 


Dffd 


[Mm] 


16.500 


Finger-Finger Spacing Thru Source 


Dffs 


[Mm] 


13.500 


Source Airbridge Inset? 


AB? 


[] 


P 


Source Airbridge Inset 


Dsabin 


[Mm] 


28.000 


Source Airbridge Height 


Hsab 


[Mm] 


3.500 


Source-Gate Airbridge Clearance 


Hgsab 


[Mm] 


1.640 


Source Pad Width 


Ws 


[Mm] 


12.000 


Drain Pad Width 


Wd 


[Mm] 


14.000 


Substrate Thickness 


Hsub 


[Mm] 


100.000 



[0053] As mentioned above, the semi-physical modeling of the intrinsic charge 

and electric field within the HEMT device is initiated by determining the relationships 
between the conduction band offset, electrical permitivities and material composition 
for the various materials in the epi stack. Material composition related band offset and 
electrical permitivity relationships may be obtained from various references, such as 
"Physics of Semiconductor Devices," by Michael Shur, Prentice Hall, Englewood 
Cliffs, New Jersey 1990. The basic electron transport characteristics, for example, for 
the linear mobility of electron carriers in the bulk GaAs cap layer may be determined 
to be 1 350cm 2 /Vs, available from "Physics of Semiconductor Devices", supra. The 
linear mobility of electron carriers in the undepleted channels is assumed to be 
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5500cm 2 /Vs. This value may be measured by Hall effect samples which have epi stacks 
grown identically to the stack in the example, except for some differences in the GaAs 
cap layer. The Schottky barrier height is assumed to be 1 .05 1 volts, which is typical of 
platinum metal on a AlGaAs material. 

[0054] The following equations represent the semi-physical analytical 

expressions to model the charge control and centroid position in the sample. 

N 5 V 



Empirical Charge Control Expression 


N, 


[cm- 2 ] 


[1+(N s 7N max )T)] 1/ r 


Ideal Charge Control with Filling Law 


Ms- 


[cm' 2 ] 


= 2N 0 ln [1+exp(V ot /(nV th ))) 


Ideal Charge Control 


No 


[cm" 2 ] 


=e. 1 V th / [2 q (dj + Adj) 1 0000] 






[cm 2 ] 


-(^maxO + NfnaxL V d8 nnmaX 


Maximum Channel Charge 




Hchan ' ^chanREF) 


Initial Gate-Channel Voltage 




tV) 


=V g .-<VAE c -V TO -aV d . 


Threshold Voltage 




m 


= O b .A Ec -V T 


Doping Threshold Voltage 




[V] 


= q N 8dalta ds 10000^! 
= {( H spac« + H bar + H fdop . + H cap )-H rec } 


Gate-to-Channel Spacing 




[m] 


(10 10 ) 




note that the expression for di can be changed for different epi-stacks 








- H chan [ 1 - d jK *V gte /H chanREF - 


Movement of Sheet Carrier Centroid 


A d i 


[m] 


diL'Vds/HchanREp] 


Empirical Charge Control Shaping Parameter 


Y 


[] 




Semi-Physical Subthreshold Populating Rate 


n 


[] 




Dielectric Permitivity of the Barrier Layer 




[F/m] 




The thermal voltage 




[VJ 


= K B T Bmb /q 


Ambient Temperature 


"l"amb 


[K] 




Fixed Em prica 1 Maximum Sheet Charge 


NmaxO 


[cm 2 ] 




vu5 ucpciiucm cinpnuai maximum oiicci oiiaiyt; 


NmaxL 


[cm 2 ] 




Vds Dependent Emprical Nmax shaping term 


n Nmax 


[] 




Channel Layer Thickness 


H chan 


[A] 




Reference Channel Layer Thickness 


^chanREF 


[A] 






(Channel Thickness for the sample for which the model was first derived) 


Schottky Barrier Height 


Ob 


[V] 




Conduction Band Offset between Channel and Barrier 




[V] 




Front Delta Doping 


Nsdeita 


[cm- 2 ] 






note that this expression can be modified for non-delta doped epi-stacks 


Gate-to-Front Delta Doping Spacing 


do 


[m] 


= {(H bar + H fdope + H^-H^} / (10 10 ) 


Barrier Thickness between front doping and channel 


u 

"space 


[A] 




Barrier Layer Thickness before front doping layer 


H bar 


[A] 




Front Doping layer thickness 


H fdope 


[A] 




Cap layer thickness 


cap 


[A] 




Empirical Drain-Induced Barrier-Lowering Term 


a 


[] 




Sheet Charge Position Gate Bias Factor 


d IK 


[A/V] 




Sheet Charge Position Drain Bias Factor 


d tL 


[AA/] 








= V th [1 + V fl /2V th + 


Effective Gate Voltage 


^gte 




sqrt(* + (V flt /2V tn -1) 2 ] 


Empirical Transition Width Parameter 


6 


[] 





[0055] As used herein, Ns represents the model sheet carrier concentration within 

the active channel. Ns' represents the ideal charge control law and is modeled as a 
semi-physical representative of the actual density of state filling rate for energy states 
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within the channel v. gate voltage. The gate-to-channel voltage used for the charge 
control, Vgt, is a function of the Schottky barrier height, conduction band offsets and 
doping in the epi stack as is known in the art. The following equations represent the 
semi-physical expression used to model the position of regional charge boundaries 
within the HEMT device. These expressions govern how to partition the model charge 
between the influence of different terminals. 





Effective Gate Length 


t-geff 




= L g + A L,+ A L d 




Gate-Source Control Region 


L .. 




= L^ + aL. + Xo, 




Source-Side Effective Gate Length Extension 


AL. 


InmJ 


= a"-.o ♦ aL k * v^ 




Drain-Side Effective Gate Length Extension 






= ALdo + ALx'V^+Al^* V^ 




Gate-Drain Control Region 


L «d 




= (LJ2 + * { tanh [10(LV2- X,,,)] +1 > / 2 


s 


















{MxdKV flt /(1+[X OL V d ,M XdL / 


■fn 


Bias Dependent Extension of the Saturated Transport Region 




Inm] 


(M XdK V gU (L ( /2 + REC 0(I ))]'"} 1 ' m 




Empirical Drain-Saturated Transport Boundary Factor 


x OL 




= L B V d ,/{2t1 + (V d A, fl . l ,ri (1/m1 } 










= X. 6 {M x .[1/(1 + {V dB /V, -n n< l ' m »- 










V^m (V^/V.^)"^ 










V.*,m[1 +(V d ,/V, -n n (1/mM> }] 




Position of the Boundary between Regions 1 and 2 


X s 




V d , M^ + V^M^} 




Note: Region 1 denotes the linear region, while Region 2 denotes the saturated region of the channel 


%sj 


Empirical Effective Gate Length Extension Gate Bias Factor 


aL k 


[ M nWJ 






Empirical Effective Gate Length Extension Drain Bias Factor 


Ak 






: ss=£ 


Effective Drain-Source Voltage Control-2 


Vdse2 


[VI 


=v d ,/[i+(v d ,/v. rtn n iym 




Rough, Intrinsic Saturation Voltage 


Vsatn 


[V] 


- It/gcM 










= 9ch 1 V at ./ 




Rough, Intrinsic Saturation Current Level 


u 


[A] 


[1 + OchiR, + sqrt(1+2g chi R J+ (V 9t ^V L f)] 




Intrinsic Conductance of the Linear Region, Under the gate 


9chi 


[S] 


={q N,H, v . Wg)/Lg 




Rough Intrinsic Saturation Voltage Level 


v t 


IV] 






Empirical Knee Shaping Parameter 


m 


[] 






Empirical Region 2 extension Drain Bias Factor 


M XdL 


[] 






Empirical Region 2 extension Gate Bias Factor 


M XdK 


[] 






Fine Intrinsic Saturation Voltage 


Vsaten 


[V] 












= gch i V flt .V t *[-V L (A+g cht R,) + 










sqa-fV^tA+g^f+V^-tg^^V,.) 2 )] / 




Fine Intrinsic Saturation Current Level 


sataom 


[A] 


rv^M i-g eW R,(v L /v Bt .) 2 )] 




Saturation Region Length Ratio 




[] 


= X 3 /L„. B 




Initial Starting position for Region 1 & 2 Boundary 






= Lg /2 




Region 1 & 2 Boundary Bias Factor 




[) 






Region 1 & 2 Boundary Drain Bias Factor 


M XsL 


[] 






Region 1 & 2 Boundary Drain Bias Factor 


M XsK 


U 





[0056] The following equations represent the semi-physical expressions used to 

model the bias dependence of linear channel mobility in depleted regions. 
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Depleted Channel Mobility 
Fixed Depleted Channel Mobility 
Depleted Channel Mobility Gate Bias Factor 



Jiave 
fidchan 
fXdK 



[cm 2 /V*s] 
[cm 2 AV*s] 
[cm 2 /V 2 *s] 



[0057] The following equations are the semi-physical expressions used to model 

the bias dependence of the saturating electric field and saturation velocity. 



Saturating Electric Field F s [V/^m] [(\i sat + R^V^) 1 0000] 

Fixed Saturating Channel Mobility (^ sat [cm 2 /V*s] 

Saturating Channel Mobility Gate Bias Factor (a satK [cm 2 /V 2 *s] 

Saturation Velocity V s [cm/s] 



[0058] FIG. 1 1 is a schematically illustrates how electrical conductance in the 

source and drain access regions are modeled in the example. 

[0059] The following equations describe the semi-physical expressions for the 
source access region conductance: 
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Source Access Resistance 



Source Access Resistance: Channel and Cap R s undepcap 



Source Access Resistance: Recess and Undeleted Cap R $ 



Col 



= (Rs, 



RsBoundaiy ) I W g 



= Rc^RF^f + 

RshI D^-tREC^+L^)] 
= Rsd. pR « : ON *MR > *tanh{[KC (K * ( 
V,. - VC^+V^MCl )) + 1 } / 2 - { V,, 

/ 2 *f1 -tanh( KR, K (V^ - VR. 0n ))]>* { tanh [ 

KR,s„ (V* - VR, Kn ..)] + 1 } / 2 + RsundepR^ * < tanhf. 

KR. K (V B ,-VR, 0n )] + 1}/2 
= Rsd.pR^° N *MR,*tanh{[KR. K *V B , + KR, L *V <ta + 
VR,ofrl + 1 } / 2 * {(1+V d ,MR, L )*MR, K *[1- 
tanhtKR.s^-VR,^.))]/ [ 2*(1 + [V g ,/((1 + 



Source Access Resistance: Crowding resistance due to conductance mismatch 


R SBoundary 




V d ,MR tL rMR, K )]^)< i n } 


Resistance of the Source Recess Access region at high on-state bias (V^) 


D ON 
"SdepRec 




= RsHd. P (REC„) 


Resistance of the Undepleted Source Recess Access region 


" S undepRec 




= RsHundap ( REC,g ) 


Uncapped, Fully Depleted Sheet Resistance 


RsHdep 


[n/sq] 


= 1 / (q N max 


Capped, Undepleted Sheet Resistance 


^SH 


fo/sq] 


= 1 / ( 1/RsHC«p + 1/ RsHundep) 


Uncapped, Undepleted Sheet Resistance 


R$Hundep 




= F .uTfUndop 1 (q N m«« *\i«deh«n) 


^ Cap Sheet Resistance 


^SHCap 


tn/sq] 


= 1 n<\ N rtap M^fRLp - H CBpEtch ) ] 


Surface Depletion Factor 


^surfUndep 


[] 




High On-state bias: Diode Turn-on voltage 


Von 




= * b -AE c -AE f 


Ohmic Contact Resistance 


R eont 






RF Ohmic Contact Resistance Reduction Factor 


R^rconF 






Source Access Resistance Bias Modification Factor 


MR, 






Cf-Vds Bias Modification Factor 


WIC^ 






Rs-Vds Bias Modification Factor 


MR $L 






Rs-Vgs Bias Modification Factor 


MR SK 






Cf-Vgs Swith point to On-state 


vc ron 






Cf-Vgs Bias Expansion Factor 








Rs-Vgs Swith point to On-state 


VRsOn 






Rs-Vgs Swith point to Off-state 


VRsotf 






Rs-Vds Swith point from Off-On transition 


VRsKnee 






Rs-Vds Bias Expansion Factor 


KR.L 






Rs-Vgs Bias Expansion Factor 


KR sK 






Rs-Vds Bias Expansion Factor @ Rs Saturation 


KR sSat 






Rs Bias Shaping Factor 


y R . 







[0060] 



The following equations describe the drain access region conductance: 
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* 



Drain Access Resistance 

Drain Access Resistance, except for recess access 
Drain Recess Access Resistance 

Drain Access Resistance: Recess and Undepleted Cap 



Resistance of the Saturated Drain Recess Access region f? osaturated 
Resistance of the Drain Recess Access region at high on-state bias (V^,) R D depRec 0N [n*n m l 



R 0 


In) 


Rdrf 


[nJ 


"dn 


tn) 


^ DAccess 




R D undapCap 




"oundepRec 





Drain Access Resistance Bias Modification Factor 
Rd-Vds Access Bias Modification Factor 
Rd-Vgs Swith point to On-state 
Rd-Vgs Swith point at saturation 
Rdaccess-Vgs Bias Expansion Factor 
Rd-Vds Saturation Bias Modification Factor 
Rd-Vds Access Resistance Bias Modification Factor 
Rd-Vgs Saturation Swith point to Off-state 
Rd-Vgs Saturation Swith point 
Rd-Vgs Saturation Swith point transition width 
Rd Bias Shaping Factor 



MR d 
MRdaccL 
VR dOn 
VRdSat 

KR daccK 

MRdSat 

MR dL 

VRdSatOfT 

VR dSa , 

: dSata 



VR, 



[] 
[] 
[] 
[] 
[] 
[] 
[] 
[] 
[] 
[] 
[] 



= ( Rorf + Ron ) 

s ( RpAcci. + RDur>d«pCap + RoundapRoc + RDSeturatad ) ' W g 

" Rdn 
= <Rdr. c )/W 9 
= R 0dopRec ON (MR d *V e . + 1)[1- 
tanh{ KF^^k • [ V,, - VR dDn + V d /MR dllceL J} ] / 2 
= R con( /RF fConF + R SH [ D gd 

- (REC Bd +L fl /2)J 
( REC gd ) 

= MR dS , t *V d , exp { -[ V 0 ,+V d /MR dL+ VR dSatOfr ] 2 / ( 2 
VRds^) } / { VR dS(lt [1+(V d ./VR dS . t ) T Rd "t"-> } 

~ RsHdep ( REC gd ) 

W26vt) + sqrt[ 5 z 

+ ((R D Ac«« " Ro«nd.pC«p)/(2SvL) - D 2 J > 



SEMI-PHYSICAL DETERMINATION OF SMALL-SIGNAL EQUIVALENT 

CIRCUITS 



[0061 ] To derive values for the familiar small signal equivalent circuit as shown 

in FIG. 1, a small signal excitation analysis must be applied to the semi -physically 
modeled physical expressions. The method of applying such an analysis is as follows: 
[0062] 1) Gate Terminal Voltage Excitation 

[0063] Apply a small+Z- voltage delta around the desired bias condition, across 

the gate-source terminals. 

[0064] Equivalent circuit element Gm = delta(Ids)/delta (Vgs') where delta 

(Vgs*) is mostly the applied voltage deltas, but also subtracting out that voltage which 
is dropped across the gate source access region, shown as RsCont, RsundepCap, 
RsundepRec, ResdepRec, and RsBoundary in FIG. 11, above. 

[0065] Equivalent circuit element Cgs and Cgd takes the form of 

delta(Nsn)/delta(Vgs)*Lgn, where delta (Nsn) is the appropriate charge control 
expression, and Lgn is the gate source or gate drain charge partitioning boundary 
length. 
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[0066] Equivalent circuit element Ri=Lgs/(Cgschannel * vs) where Cgs channel 

is the portion of gate source capacitance attributed to the channel only, and vs is the 
saturated electron velocity. 

[0067] 2) Drain Terminal Voltage Excitation 

[0068] Apply a small +/- voltage delta around te same bias condition as in 1 , but 

the delta is applied across drain source terminals. 
[0069] Equivalent circuit element Rds = l/{delta(Ids)/delta(Vds')} where Vds' 

is mostly the applied voltage deltas, but also subtracting out voltage which is dropped 
over both the gate source and gate drain access regions. 

[0070] Equivalent circuit element Cds is taken to be the sum of the appropriate 

fringing capacitance Semi-Physical models, or can take the form of 
delta(Nsd)/delta(Vds')*Xsat, were Nsd is the charge control expression for charge 
accumulation between the appropriate source and drain charge boundaries, and Xsat is 
the length of the saturated region, if in saturation. 

[0071] 3) On-mesa Parasitic Elements: The equivalent circuit elements, Rs and 

Rd are expressed by the appropriate electrical conduction models of the source and 
drain access regions. 

[0072] The RF performance can be predicted at an arbitrary bias point. 

[0073] Table 4 represents a comparison of the values for a high frequency 

equivalent circuit model derived from equivalent circuit model extraction from and 
semi-physical modeling for the sample illustrated in Table 2. 
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Table 3 

Comparison of Modeled Equivalent Circuit Results for Semi-physical 
Modeling Method, and Equivalent Circuit Model Extraction 



Hill II IdlU Utf Ul VCIld 11 


Fmiivatant Circuit 

ti|UI VCIIvlll VII villi 


Spmi-Phvsical 


Circuit Parameter 


Model 


Device Model 


Cqs 


0.227745 pF 


0.182 pF 


Rgs 


64242 fi 


infinite Q 


Cgd 


0.017019 pF 


0.020 pF 


Rgd 


133450 Q 


infinite Q 


Cds 


0.047544 pF 


0.033 pF 


Rds 


160.1791 Q 


178.1 Q 


Gm 


135.7568 mS 


124 mS 


Ri 


3.034 Q 


2.553 Q 


Tau 


0.443867 pS 


0.33 pS 



[0074] The results of the semi-physical modeling method produce a small-signal 

equivalent circuit values which are relatively more accurate than the physical device 
simulator in this case. Furthermore, given the differences in the parasitic embedding, 
treatment of the two approaches, the results given in Table 3 yield much closer results 
than a comparison of equivalent circuit values. 

[0075] Table 4 lists the values of parasitic elements used in the model 

derivations. An important difference between the extracted equivalent circuit model 
and the semi-physically derived one is the use of Cpg and Cpd to model the effect of 
launch capacitance for the tested structure. This difference leads to the results of the 
extracted model results being slightly off from the optimum physically significant 
solution. 



Table 4 

Comparison of Modeled "Parasitic" Equivalent Circuit Results for Semi- 
physical Modeling Method, and Equivalent Circuit Model Extraction 



Extrinsic Equivalent 


Equivalent 


Semi-Physical 


Circuit Parameter 


Circuit Model 


Device Model 


Rg 


1.678 Q. 


1.7 Q 


Lg 


0.029314 nH 


0.03 nH 


Rs 


1.7 a 


1.21 Q 


Ls 


0.002104 nH 


0.003 nH 
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Extrinsic Equivalent 


Equivalent 


Semi-Physical 


Circuit Parameter 


Circuit Model 


Device Model 


Rd 


3.309899 Q 


5.07 a 


Ld 


0.031671 nH 


0.02 nH 


Cpg 


0 pF 


0.02 pF 


CDd 


OpF 


0.01 pF 



[0076] As shown in FIGS. 12, 13 and 14, the modeled results that are simulated 

using the semi-physically derived equivalent circuit model very accurately replicate the 
measured high frequency, S-parameter data. 

[0077] The following equations represent the small-signal excitation derivation 

of a small-signal equivalent circuit modeled Gm. FIG. 15. illustrates the semi- 
physically simulated bias equations of the small signal Gm compared to measured data. 



Semi-Physically Modeled Drain-Source Current Control Ids 

Small-Signal Deterimination of equiv. Circuit Gm value 9 m RF 

Source-Access voltage drop ^sacc 

Fine Extrinsic Saturation Voltage Vsate 

Extrinsic Conductance of the Linear Region, Under the gate 9 C |i 
Intrinsic Conductance of the Linear Region, Under the gate 9 C hi 



[A] 
[S] 

[S] 
[S] 



=g C h v d .(inv d ,)/ 

[i + (v dB A/. ato n 1/,n 

= dl d ./d(V B .-v rtee ) 

" IdsW * ( RsUndepCap + Ac cess + RsBoundary + 
— i»atcom ^ 9cb 

[i + g C hi(Rs + Rd)] 
=(q N s ^ ave W 0 )/L 0 



[0078] The following equations represent the small-signal excitation derivation 

of Rds. FIG. 1 6 illustrates the semi-physically simulated bias-dependence of the small- 
signal Rds. 



Small-Signal Deterimination of equiv. Circuit Rds value Rds ^ 

QdsRF [S] 

Drain-Access voltage drop "dacc IV] 

Drain-Saturated Region voltage drop Vdsat W\ 

External Test probe or lead resistance RprobeD ^ 

High Frequency conductance dispersion factor r dsF [ ] 

High Frequency conductance dispersion RF r dsF [ 1 



'd,/tf(v dl 

= !<JsW * ( 



= 1 / 9d»RF 



i I 
L probeD * *ds 



> *r, 



^DUndepCap + 

R, 



'sAec" V dAce~VdSat ) 

^DUndepRec + Access + 
AAI \ 



^probeD^'g ) 
~ 'dsW * ( Saturated ) 

= (RFrd,F + D* 

tanh(10*|V - ,-V (h |)+1 
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[0079] The following equations may be used for the small-signal excitation 

derivation of Cgs and Cgd. FIG. 1 7 illustrates the semi-physically simulated bias- 
dependence of the small-signal Cgs and Cgd. 









- C gS f + (J C flcTo , Lgs* 


Small-Signal Deteri mi nation of equiv. Circuit Cgs value 




[fF^m] 


{ 1 - [ (V„ ten -V dM ) / ( 2*V sat8n -V dl J] 2 } 






= Cgdf + P C gcTot Lg d * 






[fF/^m] 


{1-[V„ ten /( 2*V„ ten -V dTO )] 2 } 


Parasitic Gate-Source Fringing Capacitance 




[fF/^m] 


= Cg surf C Mom)1 SiNF + Cgsfso urce + Cgsf Pad 


Parasitic Gate-Source Fringing Capacitance 


Cgdf 


[fF/^m] 


= Cg 8urf C Morm1 SiNF + Cgdf Cap + Cgdf Pad 


Total Specific Gate-Channel Capacitance 


CgcTot 


[fF/um 2 ] 


— Cg C + C gcdon0r 






Cg C Msh Kch(in / 


Specific Gate-Channel Capacitance 




[fF/um 2 ] 


{ [1 + (N.' / ( N 0 + J*)* 1 * 1 **]} 


Effective Drain-Source Voltage Control 


Vdse 


tv] 


=v d ./[i + (v d ,/v„ ten n i ^ 


Specific Gate-Donor Layer Accumulation Capacitance 


^ gcdonor 


[fF/ u m 2 ] 


= qdN sDonor /dV g8 


Ideal Specific Gate-Channel Capacitance 


C „c 


[fF/um 2 ] 


= qdN.'/dV flB 






[cm" 2 ] 


=<N„* + N, )*<d, + Ad,) Msh Kpar * V gte / 


Empirical Parasitic Donor Charge Control Expression 


^sDonor 


(d ( + Adj) 


Fringing capacitance to surface of source-access region 


Cgsurf 


[fF/ u m] 




Empirical Fringing capacitance-bias shaping expression 


Cf-form1 


[] 


= { 1 - tanh[ KC fK ( V g< - VC fDn + V^MC^ )] } / 


Fringing capacitance to source-access region 


gsfSourcB 


[fF/ M m] 




Fringing capacitance to source metal pads 


^ gsfPad 


[fF/ u m] 




Fringing capacitance to drain-access Capped Region 


C gdfCap 


[fF^m] 




Fringing capacitance to drain metal pads 


C gdfPad 


[fF/ u m] 




Dielectric Coating Thickness Factor 


SiNF 


[] 




Specific Gate-Channel Capacitance Bias Modification Factor 


Msh Kcahn 


[] 




Donor Charge Bias Modification Factor 


Msh Kpar 


[] 




Empirical Specific Charge Control Shaping Parameter 




[] 





[0080] The following equations are involved in the small-signal excitation 

derivation of Ri. FIG. 18, which follows, shows the semi-physically simulated bias- 
dependence of the small-signal Ri. 



Gate-Source Non-quasistatic charging resistance Rjchargo [n+n m l - L gs2 W g / [ C gsChan W g v s ] 

= C gS f + p C flC Lg S * 

Gate-Channel Source Capacitance C gsChan [fF/^m] { 1 - [ (V Mton -V dso ) / < 2*V Mton -V (to .)] 2 } 

[0081] FIG. 19 below shows the semi-physically simulated bias-dependence of 

the on-mesa parasitic access resistances, Rs and Rd. 
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EXAMPLE OF SEMI-PHYSICAL MODEL AND BIAS-DEPENDENCE 



SMALL-SIGNAL SOURCE AND DRAIN RESISTANCE. RS AND RD. 

[0082] The following example verifies how the semi-physical small-signal device 

model is able to provide accurate projections for bias-dependent small-signal 
performance. In this example, the same semi-physical device model as used in the 
previous examples was used because the example MMIC circuit was fabricated utilizing 
the same HEMT device technology. 

[0083] In this example, the bias-dependence small-signal gain and noise 

performance of a two-stage balanced K-band MMIC LNA amplifier is replicated 
through microwave circuit simulation using small signal and noise equivalent circuits 
that were generated by the semi-physical model. The results of the measured and 
modeled results are shown below in Table 5. As seen from these results, the semi- 
physical device model was able to accurately simulate the measured bias-dependent 
performance, even though the bias variation was quite wide. 

Table 5 

Measured vs. Modeled Gain NF and Gain @ 23.5 Ghz for K-band MMIC LNA at 

Difference Bias Conditions 



Bias Condition 


Measured Gain 


Predicted Gain 


Measured NF 


Predicted NF 




@ 23.5 GHz 


@ 23.5 GHz 


@ 23.5 GHz 


@ 23.5 GHz 


Vds=0.5 V 


15.2 dB 


15.8 dB 


2.97 dB 


2.77 dB 


112 mA/mm 










Vds=1.0 V 


20.6 dB 


21.0 dB 


2.29 dB 


2.20 dB 


112 mA/mm 










Vds=2.0 V 


19.8 dB 


20.2 dB 


2.25 dB 


2.15 dB 


112 mA/mm 










Vds=3.0 V 


18.9 dB 


19.1 dB 


2.30 dB 


2.11 dB 


112 mA/mm 










Vds=3.5 V 


18.4 dB 


18.5 dB 


2.34 dB 


2.18 dB 


112 mA/mm 










Vds=4.0 V 


18.0 dB 


18.0 dB 


2.37 dB 


2.27 dB 


112 mA/mm 










Vds=2.0 V 


16.4 dB 


18.0 dB 


2.45 dB 


2.21 dB 


56 mA/mm 










Vds=2.0 V 


21.4 dB 


20.9 dB 


2.38 dB 


2.21 dB 


170 mA/mm 











20 



Vds=2.0 V 


22.2 dB 


21.0 dB 


2.65 dB 


2.6 dB 


225 mA/mim 










Vds=3.0 V 


21.4 dB 


20.3 dB 


2.71 dB 


2.61 dB 


225 mA/rara 










Vds=3.0 V 


20.5 dB 


20.0 dB 


2.42 dB 


2.22 dB 


170 mA/mim 










Vds=4.0 V 


19.6 dB 


19.2 dB 


2.50 dB 


2.29 dB 


170 mA/mm 











[0084] A plot of measured vs. modeled gain for the values listed in Table 5, 
above, is shown in FIG. 20. 

[0085] The following example verifies how the semi-physical small-signal device 
model is able to provide accurate projections for physically dependent small-signal 
performance. In this example, the same semi-physical device model as used in the 
previous examples was used. 

[0086] In this example, physical process variation was input into the semi- 
physical device model in terms of statistical variation about known averages, cross- 
correlation, and standard deviations. The goal of this exercise was to replicate the 
measured DC and small-signal device variation. The degree of accurate replication 
indicates the degree to which the semi-physical model is physically accurate. 
[0087] Table 6 below lists the simulated, and known process variation that was 
used: 



Table 6 

Statistical Process Variation Model 



Parameter 


Nominal 


Standard Dev. 


Gate Length 


0.15 urn 


0.01 um 


Gate-Source Recess 


0.16 urn 


0.015 um 


Gate-Drain Recess 


0.24 urn 


0.020 um 


Etch Depth 


780 A 


25 A 


Pass. Nitride Thickness 


750 A 


25 A 


Gate-Source Spacing 


0.7 urn 


0.1 um 


Source-Drain Spacing 


1.8 um 


0.15 um 
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[0088] In the course of microelectronic component production, sample devices 

are tested in process in order to gain statistical process control monitor (PCM) data. 
FIGs. 21 A and 21 B show schematically the kind of data that is extracted and recorded 
from measured device I-V's during PCM testing. 

[0089] Since the semi-physical device model is able to simulate I-V's, it was able 

to simulate the variation of I-V's due to physical process variation. These I-V's were 
analyzed in the same fashion to extract the same parameters that are recorded for PCM 
testing. Figures 22, 23, and 24 show how accurately the simulated results match with 
measured process variation. Figure 19 shows how the semi-physically simulated Vgpk 
and Gmpk match with actual production measurements. Figure 20 shows how 
simulated Idpk and Gmpk match, also. Finally, FIGs. 21 A and 2 IB show how the 
simulated Imax and Vpo also match very well. 

[0090] Small-signal S-parameter measurements are also taken in process for 

process control monitoring. These measurements are used to extract simple equivalent 
circuit models that fit the measured S-parameters. Since the semi-physical device 
model is able to simulate these equivalent circuit models, it was able to simulate the 
variation of model parameters due to physical process variation. 
[0091] FIGs 25, and 26 show how accurately the simulated results match with 

measured/extracted process variation for the small-signal model parameters. FIG. 25 
shows how the semi-physically simulated Rds and Gm match very well with actual 
extracted model process variation 

[0092] More direct and convincing evidence supporting the accurate, physical 

nature of the semi-physical model can be shown be comparing the dependence of 
simulated and measured performance to real physical variable. As shown in FIG. 27, 
the semi-physical model is able to very accurately reproduce the dependence of Imax 
upon gate length. In addition, the semi-physical model is also able to replicate physical 
dependence for high-frequency small-signal equivalent circuits. This is shown in FIG. 
28, which shows that it is able to reproduce the dependence of Rds with Recess 
undercut width. 

[0093] Obviously, many modifications and variations of the present invention are 

possible in light of the above teachings. Thus, it is to be understood that, within the 
scope of the appended claims, the invention may be practiced otherwise than as 
specifically described above. 
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[0094] What is claimed and desired to be covered by a Letters Patent is as 

follows: 
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